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SECTION  I 


INTRODUCTION 


1 .  BACKGROUND 

— "^This  report  discusses  the  adaptation  of  a  linear  section  of  transmission 
line  to  a  time-domain  code.  The  linear  section  is  described  by  frequency- 
domain  parameters  appropriately  transformed  into  the  time  domain.  On  both 
ends  of  the  linear  section,  the  time-domain  code  employs  difference  equa¬ 
tions  to  describe  regions  which  may  exhibit  a  nonlinear  response.  The  inter¬ 
faces  between  the  linear  and  potentially  nonlinear  sections  of  cable  produce 
boundary  conditions  joining  the  variant  mathematical  descriptions  of  the 
regi  ons . — '  X 


2.  OBJECTIVES 

Shi  s  document  reports  on  the  formulation  (in  the  time  domain)  of  a 
two-port  representation  for  a  linear  transmission  line,  and  the  joining 
of  this  result  to  a  nonlinear,  time-domain,  transmission-line  code. 


In  general,  the  nonlinearities  can  occur  in  two  regions:  The  first, 
or  primary  breakdown  region,  close  to  the  burst,  results  from  the  high 
intensity  fields  associated  with  the  nuclear  device.  The  second,  or 
secondary  breakdown  region,  can  occur  at  a  discontinuity  in  the  line  far 
from  the  area  of  detonation.  These  are  caused  by  signals  propagated  by 
the  line  and  reflected  at  the  discontinuity.  These  two  regions  are  sep¬ 
arated  by  a  linear  section  of  cable. 

Of  course,  this  central  line  section  is  quite  compatible  with  the  non¬ 
linear  codes;  that  is,  the  difference  equations  representing  the  transmission 
line  equations  are  compatible  with  linear,  as  well  as  nonlinear,  line 
behavior.  .Unfortunately,  however,  the  computational  time,  and  thus  the 
computer  co6t,  is  greatly  increased,  due  to  the  large  number  of  spacial 
increments  which  must  be  devoted  to  the  linear  section.  The  question 


naturally  arises  as  to  whether  the  difference  equations  for  the  linear  sec¬ 
tion  can  be  replaced  by  a  less  time-consuming  scheme. 

As  we  shall  show,  the  linear  section  can  be  replaced  by  a  two-port 
network  section  in  the  time  domain.  The  most  obvious  formulation  leads  to 
convolution  integrals  relating  the  input  and  output  voltages  and  currents. 

The  convolution  integrals,  however,  require  an  integration  over  all  previous 
time  at  each  time  step. 

The  convolution  integrals  can  be  replaced  by  equivalent,  ordinary  dif¬ 
ferential  equations.  The  advantage  of  this  form  is  that  the  differential 
equation  solution  can  be  numerically  stepped  forward  at  each  time  step, 
based  only  on  the  value  of  the  variables  at  the  previous  time  step.  Ob¬ 
viously,  this  procedure  is  less  time  consuming  that  that  of  convolution, 
particularly  at  long  times. 

In  the  following  materials  two  approaches  will  be  shown.  The  first 
derives  the  convolution  integrals  which  relate  the  currents  and  voltages 
at  the  two  ends  of  a  linear  transmission  line,  and  the  second  approach  derives 
the  differential  equations  relating  the  same  quantities.  Because  the  dif¬ 
ferential  equations  are  most  easily  formulated  from  rational  polynomials  in 
frequency,  the  following  method  is  used  to  achieve  the  desired  representa¬ 
tion. 

The  network  parameters  are  irrational  functions  of  frequency;  so, 
using  Prony's  method,  we  approximate  the  irrational  functions  by  rational 
ones.  Then,  introducing  the  concept  of  network  state,  the  rational  fre¬ 
quency  functions  lead  to  a  set  of  first-order,  ordinary  differential 
equations  relating  the  two-port  voltages  and  currents.  These  equations 
are  in  the  desired  form. 

Finally,  the  state  equations  are  integrated  into  the  nonlinear  dif¬ 
ference  equations  to  complete  this  goal. 
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To  test  the  validity  of  the  convolution  integral  and  differential 
equation  concepts,  comparisons  are  presented  between  the  unmodified  time- 
domain  code  BLINE  and  the  two  approaches. 

In  section  III,  the  convolution  integrals  for  a  lossless,  linear  line 
are  incorporated  into  a  modified  time  domain  code  and  compared  with  the 
results  from  BLINE.  In  this  example,  field  strengths  are  limited  to  pro¬ 
duce  a  linear  result. 

The  first  example  of  the  differential  equation  procedure  is  shown  in 
section  IV.  Here,  the  transmission  line  is  a  coaxial  cable  driven  from  one 
end,  and  again,  is  contrained  to  be  linear.  Comparisons  with  the  solution 
predicted  by  BLINE  are  shown. 

The  final  example  occurs  in  section  V:  it  uses  distributed  sources 
and  the  differential -equation  approach  to  propagate  a  signal  inducing  sec¬ 
ondary  breakdown.  A  comparison  is  made  with  the  results  predicted  by  BLINE. 

All  the  examples  use  the  same  format  of  modifying  BLINE.  Specifically, 
the  linear  section  appears  as  a  single  cell  in  the  time-domain  code  inserted 
between  two  regions  described  by  difference  equations. 


SECTION  II 


THEORY 

1.  THE  TRANSMISSION  LINE  AS  A  TWO-PORT  NETWORK 

In  this  section,  we  consider  the  Laplace  transform  of  the  transmission 
line  equations  and  their  general  solution  in  terms  of  the  input  voltages 
and  currents. 

The  initial  formulation  is  in  terms  of  the  two-port,  general  circuit 
parameters  (ref.  1)  which  express  the  output  voltage  and  current  in  terms 
of  the  input  voltage  and  current.  Since,  for  this  problem,  the  admittance 
formulation  (ref.  1)  is  more  convenient,  the  equations  are  rearranged  to 
this  form.  The  admittance  formulation  expresses  the  input  and  output 
currents  in  terms  of  the  input  and  output  voltages. 

Transmission  Line  Equations 

The  Laplace  transform  of  the  single  transmission  line  equations  is 
(ref.  2) 


d_ 

dz 


Y(z,s) 

I(z,s) 


-Z(s) 1  V(z,s) 
0  I(z,s) 


e(z,s) 

i(z,s) 


(1) 


or,  in  more  compact  notation. 


y(z,s)  =  A  y(z,s)  +  g(z,s) 


(2) 


where 


0  -Z(s) 

-Y(s)  0 
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Z(s)  and  Y(s)  are  the  series  impedance  per  unit  length  and  the 
shunt  impedance  per  unit  length,  respectively, 

e(z,s)  and  i{z,s)  are  the  distributed  sources  per  unit  length, 

g(z,s)  =  [e(z,s) ,  i(z,s)], 

y(z,s)  -  [V ( 2 ,s ) ,  I(z,s)] 

s  is  the  Laplace  variable. 


The  general  solution  to  equation  (2)  is  (ref.  3) 


y(z,s)  =  e 


A(z-z J  _  f  A(z-c) 


o  ■* 


y(zQ,s)  +  f  e 

% 


g(c,s)dc 


(3) 


A(z-z  ) 

Here,  e  is  an  exponential  matrix  which  will  be  determined  shortly. 

Equation  (3)  is  the  two-port,  general  circuit- parameter  matrix  (active 
because  of  the  source  term)  for  the  transmission  line. 


Since  the  eigenvalues  of  matrix  A  are  distinct,  the  Sylvester  ex¬ 
pansion  theorem  (ref.  4)  can  be  used  to  evaluate  eA^z_z°^.  The  eigenvalues 
of  A  are  determined  by 


-A  -Z(s) 
-Y(s)  -A 


The  two  eigenvalues  of  A  are  therefore 


(4a) 


Aj  =  /Y(s)Z(s)  =  Y 
A2  =  -A(s)Z(s)  =  -y 


(4b) 
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The  Sylvester  expansion  theorm  states  that 


where 


Mz-zJ  J- 


=  £  fU,)F(A.) 
i=l  1  1 


Uz-zJ 
f(Xt)  -  e  1  0 


F(Xi)  =  j 

0=1 

j7i 


A-Xjlll 


VXj  J 


II  is  the  identity  matrix, 

Xj  are  the  eigenvalues  of  A. 

Substituting  each  equation  (4b)  into  equation  (5)  yields 


A(z-z  )  Y(z-z  ) 

e  e 

-T(z-z0) 

+  e 

1  Zo 

2 

~Yo  1 

■  a 

+  2 

Yo  \ 

y(z-zj  -r(z-zn)  -y(z-zq)  y(z-zq) 

e 


+  e 


-y(z-z0)  y(z-zq) 
e  -  e 


y(z-z  )  -y(z- 

e  +  e 


where 
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A(z-z0) 

The  elements  of  the  matrix  e  are  the  general  circuit  parameters 

for  a  two-port  representation  of  a  transmission  line  of  length  (z-zQ).  For 
the  problem  at  hand,  the  admittance  formulation  is  more  convenient,  and  can 
easily  be  obtained  by  algebraic  manipulation  of  the  variables  in  equation  (3) 
Assuming  the  line  is  of  length  i,  the  admittance  form  is 


I  (z  ,£,s)  V(£,s)  I i (z  ,t,s ) 

*  0  =  [Y]  +  1 

Iz  (zQ,Jl,s)  V(zq,s)  I2(z0,£,s) 

0  J 


where 


[Y]  = 


e-Y»._eY<'-  0 


e-yl_e-yt  o 


,-Yl  Y<  0 


yl.a-yl 
e'  +e  Y 

.-Y*  jr*  o 


I1(z0,£,s)  =  -  y0  G1(zo’*’s)  +  G2(zo’£’s) 

e  -e 


I2(z0,Jl,s)  =  _yZ_  yi  YQ  G1(zQt£,s) 
€  “6 


Gl(zo’£’s)  f°  A  (z+Z-t;) 

=  /  e  0  g(c,s)dc 

G2^zo’^’s^_  zo 


2.  TRANSFORMATION  OF  THE  ADMITTANCE  FORMULATION  TO  THE  TIME  DOMAIN 
The  admittance  formulation  of  equation  (7)  can  be  expressed  as 
IA(z0,i,s)  =  Yn(s)V(t,s)  +  Y12(s)V(zq,s)  +  I^.A.s) 

(8) 

I  (z0,t,s)  =  Y21(s)VU,s)  +  Y22(s)V(z0,s)  +  I2(z0,*,s) 
zo 

Here  Y^,  Y12,  Y21  and  Y22  are  the  elements  of  matrix  [Y] . 

Using  the  fact  that  the  product  of  the  Laplace  transforms  of  two  functions 
is  equivalent  to  their  convolution  in  the  time  domain,  equation  (7)  can  be 
directly  converted  to  the  time  domain;  that  is, 

i£(zo»A»t)  “  jf  y11(t-t,)v(A,t,)dt'  +  jf  y12(t-t,)v(z0,t,)dt‘  +  11(z0,l,t) 

(9) 

A  A 

1z(jz0»A»t^=  Jo  y21(t-t')v(£,t,)dt'  +  J  y22(t-t,)v(zo,t,)dt'  +  l2(z0,£,t) 

Here,  the  lower-case  letters  stand  for  the  inverse  Laplace  transforms  of 
their  upper-case  counterparts.  In  general,  the  inverse  transforms  of  Yll* 

Y12’  Y21*  Y22’  1V  and  *2  cannot  be  obtained  analytically,  and  must  be  ob¬ 
tained  numerically.  The  nature  of  the  inverse  transforms  for  the  admittance 
parameters  will  be  discussed  more  fully  in  section  II. 3.  The  inverse  trans¬ 
forms  of  the  equivalent  current  sources  Ij  and  I2  are  discussed  more  fully 
in  section  V. 

Equation  (9)  is  one  form  of  the  desired  current-voltage  relationship 
in  the  time  domain.  In  a  numerical  context,  however,  the  convolution  integrals 
have  a  serious  disadvantage;  because,  for  each  increment  of  time  At,  a  num¬ 
erical  Integration  must  be  done  over  all  previous  time.  Therefore,  as  time 
progresses,  the  integrations  become  longer  and  longer. 
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As  is  known  from  the  theory  of  linear  differential  equations,  a  con¬ 
volution  integral  can  represent  the  solution  of  a  differential  equation 
(ref.  3).  If  the  convolution  integrals  in  equation  (9)  can  be  replaced  by 
equivalent  differential  equations;  then,  from  a  numerical  standpoint,  the 
solution  to  the  differential  equations  can  be  stepped  forward  in  time, 
based  only  on  the  value  of  the  variable  at  that  time.  Thus,  the  computa¬ 
tional  requirements  at  each  time  step  are  greatly  reduced  from  those  re¬ 
quired  for  convolution. 

As  we  shall  see,  if  the  admittances  in  equation  (8)  are  rational  func¬ 
tions  of  frequency  (ratios  of  polynomials  in  frequency)  then  ordinary  dif¬ 
ferential  equations  can  be  obtained  which  replace  the  convolution  integrals 
(exactly)  for  all  time.  If  these  functions  are  irrational  functions  of 
frequency  (as  is  true  in  the  present  case),  differential  equations  can  be 
found  which  replace  the  convolution  integrals,  arbitrarily  closely,  for  any 
finite  time.  In  the  following  material,  we  present  a  method  to  determine 
these  differential  equations. 

Prior  to  further  consideration,  it  is  convenient  to  transfer  the  source 
terms,  1^  and  I2>  to  the  left-hand  side  of  equation  (8),  so  that 

I^(z0,«.,s)-I1(z0,J!,,s)  =  Yn(s)VU,s)  +  Y12(s)V(z0,s) 

(10) 

Iz  (zQ,fc,s)-I2(z0,Jl,s)  =  Y12(s)V(£,s)  +  Y22(s)V(zo,s) 


The  motivation  for  this  change  is  that  the  self  and  transfer  admittances 
of  equation  (10)  describe  the  response  of  the  system  initially  at  rest. 
Therefore,  the  left-hand  side  currents  are  those  which  result  from  driving 
a  system  initially  at  rest  with  voltages  V(t,s)  and  V(zQ,s). 

Consider  a  typical  term  in  equation  (10),  i.e., 


I ( s )  =  Y(s)V(s) 


(ID 
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We  assume  initially  that  Y(s)  is  a  rational  fraction  of  the  form 


liil  =  y(s)  =  ^ 
V  ( s )  D(s) 


b  sn+b„  1sn_1+...+b. 
n  n-l _ o 

a„sn+a,.  1sn"1+...+aA 
n  n-i  o 


v° 


(12) 


Note  that  the  degree  of  the  numerator  can  be  at  most  equal  to  the  degree 
of  the  denominator.  If  such  is  the  case,  then  the  impulse  response  y(t) 
corresponding  to  Y ( s )  contains  an  impulse  of  value  bn/an.  Indeed,  for  the 
transmission  line,  this  is  the  case.  It  is  possible  to  extend  this  analysis 
to  cases  where  the  degree  of  the  numerator  can  exceed  the  degree  of  the 
denominator;  however,  such  an  extension  is  not  necessary  here. 

Following  reference  (3),  we  postulate  that  the  differential  equation 
relating  i(t)  and  v(t),  in  equation  (10),  is 


Jr  «t>  *  Vi 


(13) 


v(t)  +  b 


n-l 


d"'1 

dt"'1 


v(t)  +*--+bQv(t) 


The  Laplace  transform  of  equation  (13)  should  be  reducable  to  equation  (12). 
The  Laplace  transform  of  equation  (13)  is 

an  {snI(s)  -  sn"11(o)  -  sn'2  ^  i(o) . ^  i(o) 

d  t 

♦  Vl  K'1'!5)  -  -  s"-3  3t  '(")  --  ‘^2  «»>} 

+  •  I  • 
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=  bn|sn  V(s)-sn'1v(o)-sn‘24_v(o) - ~-y  v(o)| 

+  bn-l  {sn'1V(s)-sn'2v(o)-sn-3^v(o)-...  ^2  v(o) 


Solving  equation  (14)  for  I(s),  the  result  is 

I(s>  *  pfff  V(S)  *  ^  |sn_1  [ani(o)-bnv(o)] 

*  s"'2  [an  3t  ,(o)  -  bn  at  v(o>  *  an-li(o)  ‘  bn-lv(o>] 

(15) 

*[an  dt571” 1<0)  '  b"  at"71  V<0>  +-',+ai1(o)-biv(o)]} 
Here  N(s)  and  D(s)  are  defined  in  equation  (12). 


If  the  transfer  function  of  the  system  is  to  be  that  of  equation  (12); 


i  .e. , 

i<s)  *  Irff}  v<5> 

then  all  other  terms  in  equation  (15)  must  be  equal  to  zero.  These  con¬ 
ditions  form  the  interrelationship  among  the  initial  conditions  on  the 
variables.  However,  rather  than  using  these  initial  conditions  directly 
(v(t)  is  the  source  and  i(t)  the  response),  it  is  more  convenient  to  make 
use  of  the  concept  of  system  state. 


■  zji, 


Basically,  the  state  of  a  system  is  a  property  of  the  system  which, 
together  with  system  inputs,  determine  the  future  state  of  the  system. 
For  present  purposes,  the  state  can  be  a  set  of  new  variables  defined  as 
linear  combinations  of  the  original  variables,  i(t)  and  v(t),  and  their 
derivatives.  Since  the  transmission  line  is  initially  at  rest,  and  the 
fact  that  we  wish  equation  (15)  to  follow  from  equation  (14),  a  very  con¬ 
venient  definition*  of  the  state  variables  x^t)  is 

xx(t)  =  ani(t)-bnv(t) 

x2(t>  ■  an  at  1<t)_bn  3t  (16) 


x„(t>  ■  an  ~iPf  iU)'bn  £^V(t)  t"-+al1<t>  -  blv(t> 


It  Is  easily  seen  that  if  the  initial  state  (xj(o),  x2(o) ,. . . ,xn(o))  is 
zero,  then  equation  (11)  will  result  from  equation  (15).  We  have  thus 
arrived  at  a  set  of  variables  whose  initial  value  is  zero,  and  which  re¬ 
sults  in  the  desired  relationship. 

We  now  determine  the  differential  equations  which  these  new  state 
variables  must  satisfy. 

Equation  (16)  can  be  written 


xx(t)  -  ani(t)  -  b„v(t) 

-2(t)  *  at  <!<*)  *  vi«t)  -  Vi,(t>  on 


x„(t) 


d_ 

dt 


xn_i(t)  +  ax1(t)  -  bjV(t) 


*  System  state  is  not  unique. 
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The  last  state  variable  xn(t)  in  equation  (16)  can  be  combined  with 
equation  (13)  to  yield 

3t  xn(t)  =  ’  {aoi(t)  •  bov(t)}  08) 

Finally,  solving  equations  (17)  and  (18)  in  terms  of  ^  xfc(t)  and 
i(t),  and  using  the  first  relationship  in  equation  (17);  i.e., 

i  bn 

i(t)  =  —  x,(t)  +  v(t), 

n  n 

leads  to  the  result 

H  xl(t)  {-a„-ix1(t)*arix2(t)»(a„b„.i-a„.1b„)v(t)} 

3t  x2(t)  ’  ^  {-an-2xl(tl+V3(t)+(anbn-2-V2t'n>v<t>) 

•  (19) 

3t  Vl(t)  =  falxl(t>tanxn(t)+(anbrall’r,)v(t)} 

at  xn(t>  =  5;{-aoxl(t)+(anbo-aobn)v(t)} 
i(t)  =  a-  {xi(t)+bnv(t)} 

Equation  (19)  is  a  set  of  first-order,  ordinary,  differential  equations, 
which  with  zero  initial  conditions,  describe  the  desired  voltage  current 
relationship  between  i(t)  and  v(t). 
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3.  IRRATIONAL  ADMITTANCE  ELEMENTS 


In  the  previous  section,  a  technique  is  described  which  replaced  a 
convolution  integral  with  a  set  of  ordinary  differential  equations.  The 
method  is  useful,  however,  only  if  the  individual  admittance  elements  are 
rational  functions  of  frequency. 

For  the  transmission  line,  the  admittance  elements  are  irrational 
functions  of  frequency,  so  it  is  necessary  to  approximate  them  by  rational 
ones,  if  this  method  is  to  be  used.  There  are  many  ways  to  make  this 
approximation  (see,  for  example,  reference  5).  A  particularly  convenient 
method  and  the  one  that  is  used  here,  is  Prony's  method  (ref.  6). 

Using  Prony's  algorithm,  it  is  possible  to  represent  a  time  function, 
arbitrarily  closely  over  any  finite  time,  by  a  finite  sum  of  decaying 
sinusoids*.  Thus,  the  representation 
N  $  t 

f(t)  A.  e  1  (20) 

i=l  1 


is  obtained.  The  Laplace  transform  of  this  expression  is 


F(s)  =  =  E  A-  .  "U 


s-s  • 

1=1  1 


which  is  a  rational  function. 

To  use  the  Prony  procedure  it  is  necessary  to  inverse  transform  the 
admittance  functions.  For  the  transmission  line,  most  of  the  admittance 
functions  are  sufficiently  complex  that  the  inversion  must  be  done  numer 
ically.  In  addition,  these  Inverse  transforms  contain  Impulses.  The 
value  of  these  impulses  must  be  removed  before  numerical  inversion. 

*  Simple  decaying  exponentials  are  included. 


To  illustrate  the  procedure,  consider  the  characteristic  admittance 
of  a  transmission  line  with  lumped  series  resistance  and  inductance  and 
only  shunt  capacitance  (perfect  dielectric).  The  characteristic  admittance 
for  this  line  is 

Y  /o  -  JUT  (22) 

Vs'  IR+sL 
We  note  that 

lim  Y(s)  = 

and  since  the  Laplace  transform  of  an  impulse  is  a  constant,  the  inverse 
transform  of  Yq( s )  contains  an  impulse  of  value  ^C~ .  Therefore,  if 

we  were  to  numerically  inverse  transform  YQ(s),  we  would  remove  this  im¬ 
pulse  and  inverse  transform 

y;<s>-y„<*> -VF-Vk^-VF  U3) 

The  prony  procedure  would  then  be  applied  to  the  result.  Of  course,  an 
impulse  of  value  /(T  would  then  be  added  to  the  exponential  approximation 

Ml 

of  equation  (20). 

It  turns  out  that  equation  (22)  can  be  analytically  inverse  trans¬ 
formed,  however,  and  the  result  is 


V2L 


(24) 


The  impulse  occurs  in  equation  (24),  because  of  the  unit  discontinuity  of 
the  expression  in  parenthesis  at  the  origin. 
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Equation  (24)  cannot  be  entirely  represented  by  a  sum  of  simple  expo¬ 
nentials.  However,  a  Prony  procedure  applied  to  the  first  30  microseconds 
of  this  function  With  the  values 

_2 

R  =  1.3611  x  10  ohms/meter 
L  =  5.8322  x  10‘8  henry/meter 
C  =  4.3878  x  10"10  farad/meter 

indicates  that  two  simple  poles  form  a  good  approximation. 

The  following  poles  and  residues  result  from  the  Prony  analysis: 

Pole  Residue 

-5.7362xl04  1.826xl03 

-2.0056xl05  8.298xl03 

A  comparison  of  the  resulting  approximate  function  with  the  original  ir¬ 
rational  function  is  shown  in  figure  1.  Actually  all  values  are  within 
.1%  except  for  the  last  several  values  which  are  within  1%. 

The  approximation  to  YQ(t)  is  then 

Y0(t)  =  <5(t)+l . 826xl03Exp(-5. 7362xl04t)+8.298xl03Exp(-2.0056xl0St) 

Here  the  impulse  has  been  added  to  complete  the  expression.  The  Laplace 
transform  of  this  function  will  form  a  rational  approximation  to  Yq( s ) . 


NORMALIZED  AMPLITUDE 


I 


TIME  (ps) 


Figure  1.  Comparison  of  Irrational  Function  With 
A  Two  Exponential  Approximation 
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4.  ANALYSIS  OF  THE  ADMITTANCE  MATRIX 


It  is  instructive  to  analyze  the  admittance  matrix  [Y]  in  equation  (7) 
in  more  detail,  by  expanding  the  individual  elements  in  the  following  way: 


eYVe"Y£  Y  =  _  1+e"2^ 


e-y£_ey£  o  i_e-2Y£  o 

=  -(l+e"2Y£)(l+e"2Y£+e'4YV--)  Y£ 
=  -(l+2e'2Y£+2e“4YS,+  ---)  Y Q 

■Tt2  yn  \  s  -Ze-^d+e-^^e-4^...)  Y 
e  -e 

=  -2(e“Y£+e“3Y£+e“5Y£+*--)  Y 

o 

Therefore,  the  admittance  matrix  can  be  expressed  as 


(25a) 


(25b) 


[Y]  = 

-1  O' 

Y  +2e"y£ 

0 

■  o  r 

Y  +2e”2y£ 

0 

-1  o ' 

,  0  1. 

.-1  0. 

.  0  1  . 

V 


(25c) 


We  note  that  the  multiplicative  factors  of  2  in  this  expression  are  due  to 
the  fact  that  this  is  a  short-circuit  admittance  matrix,  and  the  current 
reflection  coefficient  under  short-circuit  conditions  is  2. 

Due  to  the  multiplicative  factors  of  e~nY£(n=0,l,2,. . .)  in  equation  (25), 
each  matrix  in  the  sum  undergoes  a  delay  (in  the  time  domain)  of  multiples 
of  a  line  length  delay.  If  the  line  is  lossy,  this  delay  is  of  course 
dispersive;  nevertheless,  it  is  there.  Therefore,  the  problem  can  be  sub¬ 
divided  (time-wise)  into  a  number  of  distinct  parts,  each  of  which  is 
additive  to  the  total  result.  Furthermore,  since  we  are  usually  interested 


22 


'  I  ■•*■1  <(****&*■** 


Y  =\s  LC  =  s  VtC 


so  that  the  exponential  factors  represent  only  time  delay;  i.e.,  multi- 

-st 

plication  of  a  Laplace  transform  by  e  0  delays  the  ensuing  time  function 

by  V 

Since  Yq  =  is  constant,  its  inverse  transform  is  an  impulse  of 
value .  Therefore,  the  inverse  transform  of  equation  (25)  is 


-6(t-2t0) 


0 

6(t-2t 


k* 


(26) 


In  this  equation,  t  =VTC  l  . 

Each  term  in  equation  (26)  is  a  simple  impulse,  so  that  the  convolution 
integrals  in  equation  (9)  are  easily  evaluated.  The  result  is 

iU»t)  .VF{-vU,t)«v(z0.t-t0)-2v(£.t-2tb)  +  }  +  i x ( «- , t ) 


i(z0,t)  =^"{v(z0,t)-2v(S,,t-t0)+2v(z0,t-2t0)  -  •••}  +  i2(£,t) 


(27) 


Equation  (27),  then,  is  the  current-vol tage  relationship  for  a  lossless 
line.  In  this  expression  all  voltages  are  zero  if  t-kt  is  negative.  In 
section  III  a  different  derivation  of  the  voltage  and  current  relationship 
and  numerical  results  are  presented  for  a  lossless  transmission  line. 

When  the  line  is  not  lossless,  the  inverse  transform  of  equation  (25c) 
is  more  complex.  In  general,  each  term  contains  an  impulse  plus  a  dis¬ 
persive  term,  so  that  the  convolution  integrals  must  be  numerically  eval¬ 
uated. 

5.  BURIED  CABLE  CONSIDERATIONS 

The  motivating  problem  for  this  study  requires  that  we  determine  the 
voltage  across  the  dielectric  of  a  buried,  insulated  conductor;  that  is, 
we  must  determine  the  dielectric  voltage  to  determine  dielectric  breakdown. 
In  the  previous  sections,  we  have  analyzed  the  behavior  of  a  single  trans¬ 
mission  line,  with  no  special  consideration  for  the  dielectric  portion  of 
the  total  line  voltage.  It  is  the  purpose  of  this  section  to  present  the 
procedures  to  be  followed  to  determine  the  dielectric  voltage. 

It  is  common  practice  in  the  analysis  of  buried  cables  to  represent 
the  dielectric  and  ground  in  the  way  illustrated  in  figure  2.  Typically, 
and  Yg  each  consist  of  a  parallel  resistor-capacitor  combination;  so 
that,  in  the  time  domain,  the  transmission  line  equations  can  be  written 


a(vd+vc 

3Z 


-  Ri  +  E1nc(z,t) 


3i  _  r  r  „ 

3z  "  ”  cd  3t  "  bdvd 


3i 

3z 


3v„ 

C  — 2.  _ 
3t 


G9V9 


(28) 
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Here,  the  subscripts  d  and  g  refer  to  dielectric  and  soil,  respectively. 


UNIT  LENGTH  OF  BURIED 
TRANSMISSION  LINE 


Z  =  SERIES  IMPEDANCE/UNIT  LENGTH 

Yd  =  SHUNT  ADMITTANCE/UNIT  LENGTH  OF  THE  DIELECTRIC 

Y  =  SHUNT  ADMITTANCE/UNIT  LENGTH  OF  THE  GROUND 


Figure  2.  Section  of  Buried,  Insulated  Conductor 

Indeed,  the  present  nonlinear,  time-domain,  difference  equations  are  re¬ 
presentations  of  these  partial  differential  equations. 

To  proceed,  we  again  resort  to  the  Laplace  transform  versions  of  the 
single  transmission  line  equations  which  are 


|L  =  .  yv 

az 


Here,  however,  to  conform  to  figure  2,  the  line  parameter  Y  is  the  series 
combination  of  Yd  and  Y^;  i.e., 

1 

Y 

O  Y 

Now  the  dielectric  voltage  can  be  written  in  terms  of  line  voltage 
as  follows: 

1 


Therefore,  in  terms  of  dielectric  voltage  Vd>  the  total  line  voltage  V  is 
Y  +  Y 

V  =  d-  y -  3  vd  (30) 

g  1 

Similarly,  the  line  voltage  V  can  be  expressed  in  terms  of  ground  voltage 
Vg;  that  is, 

Y,  +  Y 

V  >  »  (31) 

d  y 

Expressions  (29)  and  (30)  are  the  desired  relationships  to  convert  from 
Vd  or  Vg  to  V  and  V1'ce-Versa-  However,  since  V  =  Vd  +  Vg,  it  is  only 
necessary  to  use  one  of  these  expressions.  Both  of  these  expressions  are 
Laplace  transforms,  so  that,  in  general,  the  methods  of  parts  3  and  4  of 
section  must  be  used  to  obtain  the  corresponding  time  domain  relationships. 

We  mention  in  passing  that  the  transmission  line  equations  can  be 
written  solely  in  terms  of  dielectric  voltage  Vd.  This  is  accomplished 
by  substituting  equation  (29)  into  the  transmission  line  equation  (28); 
the  result  is 


I_  +  i 

d  g 


=  IA 

YW+Yn 

d  g 
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r 

££ 


h 

& 


V 

V 

f 

I 


l 


!!d 

az 


V  ■  +  Y 
d  g 


-Z'l 


31 

az 


-Ydvd 


Y  1 

where  V  =  y — S__ 
d  g 


(32) 


Equation  (31),  while  having  different  line  parameters  than  equation 
(28),  has  the  same  propagation  constant.  Therefore,  the  linear  section 
of  transmission  line  can  be  modeled  either  with  total  voltage,  or  with 
dielectric  voltage. 
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SECTION  III 


LOSSLESS  LINEAR  SECTION 


1.  LIMITATIONS  OF  THE  APPLICATION 


This  section  of  the  report  uses  the  straightforward  approach  of  con¬ 
volution  integrals  (equation  9)  to  find  the  relationship  between  the 
currents  and  voltages  at  the  ends  of  the  linear  insert.  With  no  distri¬ 
buted  sources,  equation  (9)  becomes 


/ 


i(z„,t)  =  /  yu(t-f)vU.t 


y12(t-t')v(zo,t)dt'  (33a) 


t  c 

i U,t)  =  f  y21(t-t')v(«,,t ')dt'+  J  y22(t-t’)v(zo,t)dt ' 


(33b) 


To  calculate  the  admittance  matrix  as  a  function  of  time,  the  inverse 
Fourier  transform  of  the  frequency  representation  must  be  found.  As  dis¬ 
cussed  in  section  II. 3,  both  the  self  and  transfer  admittances  contain  im¬ 
pulses;  as  a  result,  if  the  impulses  are  not  removed,  the  inverse  Fourier 
transform  is  not  defined.  Another  method  to  find  the  current-voltage 
relationship  is  the  following.  Examining  only  the  first  integral  in 
equation  (33a), 


yn(t-t')vU,r)dt 


F_1(yi:l(s)V(£,s)) 


-  F_1[(“T~)  (sV(£,S))]=  F_1  Yu(s)V(£,s) 


(34) 
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The  quantity  Y^s),  has  no  impulse  behavior  and  therefore  a  well  defined 
inverse  Fourier  transform;  it  contains  instead,  discontinuities  in  the  time 
domain.  The  voltage  part  of  the  integral,  V(i,,s),  simply  represents  the 
time-domain  derivative  of  the  voltage.  The  equation  now  appears  as 

Jt  t  _ 

J  yn(t-t ')vU,t ')dt'  =  f  yu(t-t')  v(£,t')dt '  (35) 

O  (0 

The  other  elements  of  the  admittance  matrix,  y12(t),  y21(t),  and  y22(t) 
will  have  properties  similar  to  y^(t);  the  following  discussion  of  y^t) 
can  be  applied  to  all  other  admittance  terms  and  their  contributions  to  the 
currents. 

The  time  behavior  of  y^( t)  consists  of  a  discontinuity  at  t  =  0 
(t  =  Yl,  one  delay  length  for  the  transfer  admittance  terms)  and  at  each 
subsequent  2nYJL  Between  the  discontinuities  at  the  reflection  times, 
y^(t)  is  monotonically  increasing.  If  the  line  is  lossy,  then  in  the  limit 
as  t-*<»,  y^t)  =  constant.  If  the  line  is  not  lossy,  the  monotonic  sections 
are  flat,  the  discontinuities  do  not  diminish  with  time,  and  there  is  no 
limit  to  y^t). 

The  only  case  analyzed  with  the  admittance  matrix  in  this  form  was 
for  a  lossless  line;  the  size  of  the  discontinuity  was  then  a  constant 
equal  to  the  characteristic  admittance  of  the  line  (occuring  at  t  =  0  and 
each  two  line  delay  lengths).  In  this  form  it  was  possible  to  Include  all 
reflection  terms  and  signals  propagating  in  both  directions. 

2.  ADAPTATION  TO  BLINE 

The  junction  of  the  linear  section  with  the  sections  described  by  dif¬ 
ference  equations  employs  the  solution  of  the  dielectric  and  ground  voltages 
at  the  junction  points.  The  transmission-line  equations  (equations  28) 
employed  by  BLINE  are 
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i; 


3z^vd+vg^  =  "L  3t  _Rc1 

(35a) 

3i_  f  3Vd  rv 

3z  "  "Cd  at  ■  Gdvd 

(35b) 

vr  =  -C  Kf"  -  G  V 

3z  g  3t  g  g 

(35c) 

when  no  sources  are  included. 

An  examination  of  figure  3  shows  the  spatial  difference  equation  for 
the  current  employing  the  self  and  transfer  impedances  is 

M.  =  iL_  ( r  _  j  ) 

3z  Az  uj+2  ‘j+l+V 

y21(t-t‘)  dt>  (36a) 

3v(z  .f)  \ 

y22(t-t,}  —JV -  dt  ) 

on  the  right  end  of  the  linear  section  (j+2)  and 


iLL  =  iL  (j  _  i  ) 
3z  AZ  Uj+ls  ly 


3v(z  .f) 


.a  f-y  (t-t'l  flifal 

az  \ y  yir  z  ‘  dt  • 


dt' 


A 


(36b) 


30 


'1  'l  v2  I; 


-I  *  I 


Linear  Insert 


V4  I ,  v 


j  Aj  >1 


Vj+2  1  j+2  Vj+3 


Figure  3.  Placement  of  Current  and  Voltage  Cells  in 
BLINE  and  the  Linear  Insert  in  BLINE 


on  the  left  end  (j+1).  In  both  cases  v  is  the  total  voltage  as  given 
in  section  II. 1. 

Approximating  the  convolution  integrals  as 


3v(z0,t') 

9p 


k 


(37) 


where  the  k  and  n  indices  represent  time;  equations  (36a),  (35b),  and  (35c) 
then  give  the  boundary  conditions  on  the  right  side  of  the  linear  section 
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where  the  quantities  subscripted  by  d  and  g  indicate  capacitances,  conduc¬ 
tances  and  voltages  across  the  dielectric  and  ground  and 


where 


The  resulting  boundary  conditions  are  simultaneous  equations 
and  vj*> 

3.  COMPARISON  WITH  BLINE 


(39c) 


in  V 


n+1 

d 


The  comparison  problem  run  by  BLINE  and  BLINE  modified  by  the  convol¬ 
ution-integral  approach  consisted  of  a  201  meter  cable;  the  linear  insert 
consisted  of  the  center  one  meter.  The  other  parameters  are  shown  in 

in  table  1. 


TABLE  1 


rl 

r2 


rd 


Core  radius 

Dielectric  radius 

Soil  dielectric  constant 

Soil  conductivity 

Sheath  dielectric  constant 


2.5xl0~2m 

2.627xl0"2m 

10  farads/meter 

_3 

10  mhos/meter 
2.3  farads/meter 


The  observation  points  used  in  both  analyses  are  at  the  two  ends  of 
the  cable.  The  line  was  pulsed  by  a  voltage  excitation  at  the  left  end 
of  the  line  with 


Figure  4  shows  the  total  voltage  as  a  function  of  time  predicted  by 
BLINE  and  BLINE  with  the  modification  of  the  1  meter  linear  section. 

Figure  5  compares  the  total  voltage  at  the  end  of  the  cable  opposite 
from  the  source. 

The  voltage  comparison  at  the  source  end  of  the  cable  shows  that  no 
reflection  occurs  at  the  interface  between  the  linear  and  difference- 
equation  described  sections.  The  voltages  at  the  far  end  of  the  cable 
show  a  lower  amplitude  after  the  pulse  travels  through  the  linear  section 
when  compared  to  the  unmodified  BLINE  predictions;  the  cause  of  this  dis¬ 
crepancy  is  presently  unknown. 
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SECTION  IV 


LINEAR  SECTION  IN  A  COAXIAL  LINE 

1.  LIMITATIONS  OF  THE  APPLICATION 

The  reason  for  the  insertion  of  a  linear  section  into  the  analysis 
is  to  save  the  computer  time  necessary  to  propagate  the  signal  in  the 
linear  section.  In  this  region,  the  time  domain  analysis  is  not  necessary 
since  no  breakdowns  are  expected;  parameters  linking  the  voltages  and 
currents  at  the  ends  of  the  linear  section  are  ail  that  are  necessary. 

Also,  since  the  prediction  of  secondary  breakdown  phenomena  is  desired, 
one  actually  needs  only  the  signal  propagating  from  the  source  through 
the  linear  section.  In  the  region  of  secondary  breakdown,  this  signal  is 
used  to  drive  the  nonlinear  section  and,  analyzed  by  the  time  domain  code, 
can  reflect  from  the  end  of  the  line  causing  nonlinearities.  This  reflected 
signal,  returning  to  the  source  region,  cannot  reasonably  be  expected  to 
traverse  the  entire  cable  and  add  nonlinearities  to  the  primary  break¬ 
down  region;  in  this  analysis  its  presence  in  the  linear  region  is  ig¬ 
nored. 

Examining  only  the  first  pulse  from  the  source  end  of  the  line,  the 
current  and  voltage  relationships  from  equation  (25c)  are 

IU.S)  =  -Yq  VU.s)  +  2e‘n  Yq  v(zQ,s)  (41a) 

Kz0,s)  =  Y0  v(zo»s)  (41b) 

This  form  is  used  since  all  subsequent  admittance  terms  in  equation  (25c) 
represent  delay  lengths  associated  with  signals  propagating  toward  the 
source  and  their  reflections.  The  terms  of  the  admittance  matrix  which 
are  present  are:  the  self  admittance  for  both  I(£,s)  and  I(z0,s)  and  the 
transfer  admittance  delayed  by  one  line  length  for  I(£,s). 
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The  analysis  of  a  coaxial  line  was  chosen  as  the  next  step  in  order  to 
introduce  losses  into  the  linear  section  via  a  nonzero  resistance  of  the 
center  conductor.  The  shield  was  chosen  to  have  a  zero  resistance  so  the 
entire  system  would  be  represented  by  a  single  voltage;  ordinarily,  the 
calculation  of  two  voltages  are  necessary,  one  across  the  dielectric  and 
one  across  the  soil.  Inserting  a  large  soil  conductivity  into  the  time 
domain  code,  resulted  in  the  approximation  of  a  coaxial  cable  with  no 
shield  resistance  and  no  voltage  drop  across  the  soil. 

Constant  R,  L,  G,  and  C  (section  1.3)  parameters  were  also  used  in 
the  analysis  to  allow  the  subsequent  steps  in  the  development  of  the  pro¬ 
cedure  to  be  more  easily  understood. 

2;  PRONY  REPRESENTATION  OF  THE  ADMITTANCES 

The  representation  of  Yq  by  an  impulse  term  and  decaying  sinusoids 
was  shown  in  section  1.3.  The  procedure  to  find  the  representation  of  the 
transfer  admittance  term  is  identical.  First,  the  impulse  must  be  identified 
and  handled  appropriately. 

Assuming  the  coaxial  cable  contains  only  a  shunt  capacitance  term,  the 
characteristic  admittance  and  propagation  constant  are 


vs>=VrS?  (42a) 

Y(s)  =VsCTR+sLT  (42b) 

At  high  frequencies,  the  limits  of  these  are 
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Tim  Y  (s)  =Jf  (43a) 

S-K»  1 

lim  Y (s)  =  l  R  Jp  -  sVLC1  (43b) 

Thus  the  impulses  of  the  admittances  are 

Impulse  (Yq)  = 
and 

Impulse  (2e"y£Yo)  =  exp  [  -  (|  R^p-  s  VLC1)  & 

Next,  the  Prony  process  is  applied  to 

f_1[yo(s)-VF]  (45a) 

and  the  transfer  admittance  advanced  by  one  delay  length, 

F’1  j  |2e~ya  Yq  -  2^p  exp(-{lRVu"  s  VlT} n)]e‘s [  (45b) 

It  is  worthwhile  to  notice  here  that  the  Prony  results  of  the  self 
admittance  term  are  independent  of  the  line  length,  whereas  the  transfer 
admittance  term  is  a  function  of  the  line  length.  The  Prony  analysis  is 
done  on  the  transfer  admittance  advanced  by  one  delay  length  to  reduce 
the  number  of  poles  necessary  to  achieve  an  adequate  fit.  The  inverse 
transform  of  the  transfer  admittance  is  near  zero  until  one  delay  length 
in  time;  at  that  time,  the  admittance  is  discontinuous  to  a  finite  value 
after  which  it  decays.  This  type  of  time  function  can  be  represented  by 
a  sum  of  decaying  exponentials,  but  far  fewer  poles  are  required  if  the 
discontinuity  in  the  funciton  is  shifted  back  to  the  time  origin. 


(44a) 

(44b) 
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Another  difficulty  arises  in  the  inverse  transform  of  a  discontinuous 
function;  Gibbs'  phenomenon  is  present  in  the  transform  resulting  in  an 
inaccurate  answer  at  the  time  origin.  Two  approaches  are  used  to  reduce 
this  effect:  the  first  is  to  calculate  the  frequency  domain  function  to 
frequencies  far  higher  than  would  appear  necessary  to  resolve  the  actual 
behavior,  f^lO11  Hz.  This  increases  the  Gibbs'  oscillation  frequency  and 
decreases  the  decay  time.  The  second  method  is  to  analytically  calculate 
the  inverse  transform  at  t  =  0.  This  can  be  accomplished  by  taking  the 
limit, 

-  f(o+)  (46) 

the  inverse  transform  of  F(s).  These  limits  for  the  ad- 

[ s  (-Vs> +  VF)]-  fr  VF  <47a> 

js  f-2e'Y\(s)  +  2 ^  exp  [-(£  -  sVLC1  l\]\ 

1  1  (47b) 

exP(svt? 

Again,  the  delay  time.VLC  l  shows  in  the  limit  of  the  transfer  admittance. 

Using  the  limit  as  s  goes  to  infinity  of  sF(s)  as  the  value  of  f(o), 
the  two  admittances  were  approximated  by  the  Prony  program  as 

F-1  [  VF'  Y„]=  a/’'  +  ^  <48a> 


lim 

S-M» 


sF(s) 


where  f(t)  is 
mittances  are 


lim 

$-x» 


and 


1  im 

S-*» 
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or 


vs)  *  VF-  5^ 

a2 

S-S2 

(48b) 

where 

ax  =  1.82602x10^ 

a2  = 

8.29867xl03 

Sj  =  -5.7362xl04 

s2  = 

-2.00559x10^ 

and 

F"1  j-2e_Y£Y0(s)  + 

2^p  exp  [  -  | 

(iR^1  -  sVcrKIj 

1  sl* 

*aj  e  +  a. 

s2t 

2 

(49a) 

or 

2e'Y£  Yq(s)  =  2^ 

exp  [-(  j  Ri 

al 

a2 

(49b) 

s-sl 

s-s2 

where 

aj  =  1.9269xl03 

a2  = 

5. 98615xl03 

Sj  ■  -4.68392xl04 

S2  = 

-1.93015X105  • 

The  values  of  the  two  impulse  terms  were 

Impulse  (Yq)  =  8.6737xl0"2 
and 

Impulse  [  2e'y£  Y  (s)]  =  9.6134xl0-2 

In  both  cases,  a  fit  using  two  poles  and  residues  achieved  an  error  of 
less  than  one  percent  at  all  except  the  last  two  time  points. 


41 


The  admittances  can  now  be  represented  by  ratios  of  polynomials  in  s 


v  -  Kll5  *  W12s  *  "13 
Dlls2  +  °125  +  D13 


and 


2e'Y\(s)  = 


N21s  +  N22S  4  N23 
D21s  +  °22S  +  °23 


(50a) 


(50b) 


where  the  numerator  and  denominator  in  both  cases  are  of' the  same  order  in 
s  because  of  the  impulses  present  (section  II. 2).  The  transfer  admittance 
is  understood  to  be  delayed  in  time  by  one  line  length. 

With  the  ratios  of  polynomials  known  for  the  two  admittances,  the  state 
equations  can  be  derived.  From  section  II. 2,  the  state  equations  and 
initial  conditions  are 


xu(0)  =  0,  x12(0)  =  0 

x2 i ( 0 )  =  0,  *22 ( 0)  =  ^ 

•  ^12  /  °12  \ 
xn  =  ■  X11  +  X12  +  \  12  "  N11  DJY/  V 

•  ^13  /  ^13  \ 

x12  =  '  X11  +  v  13  '  N11  )  V 

x2i  =  "  0^  X21  +  x22  +  (N22  ‘  N21  )v 

•  /  ^23  1 

x22  =  "  X21  +  \  23  ‘  N21  D”  )V 


(60a) 

(60b) 

(60c) 

(60d) 
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where  the  first  subscript  refers  to  the  self  admittance  (1)  or  the  transfer 
admittance  (2).  The  current  as  a  function  of  the  voltage  is  now  given  by 


*11  Nii 

i  =  =-=-  +  = —  v.  (self  admittance)  (61a) 

D11  D11 

for  the  self  admittance  part  and 

*?1  Npi  _ 

i  =  jr— =-  +  n —  v  (transfer  admittance)  (61b) 

°21  U21 

for  the  transfer  admittance  portion.  The  quantity  v  represents  the  voltage 
delayed  by  one  propagation  length  of  the  linear  section.  Because  there  is 
interest  in  only  the  signal  propagating  from  the  source  in  the  linear  portion, 
the  voltage-current  relationship  is  given  by  only  eauation  (61a)  for  the  end 
of  the  linear  section  nearest  the  source.  At  the  other  end  of  the  linear 
section,  the  voltage-current  relationship  is  represented  by  the  sum  of  the 
two  current  contributions. 

3.  ADAPTATION  TO  BLINE 

The  reason  for  the  formulation  of  the  equations  giving  current  as  a 
function  of  voltage  was  compatibility  with  BLINE.  The  boundary  conditions 
in  BLINE  always  employ  a  current  as  a  function  of  voltage. 

The  solution  of  the  state  equations  is  easily  realizable  in  BLINE 
using  the  formulation 

x( t+At)  =  x(t)  +  x(t)At.  (62a) 


This  gives 


.  Oi  o  . 

Xn(t+At)  =  Xn(t)  ^  1  -  At  j+  X12(t)  At 


r\ 


(62b) 


(62c) 


for  the  transfer  admittance  current  contribution. 

In  the  case  of  a  coaxial  cable,  the  voltage  across  the  shield  is  zero 
and  only  the  dielectric  voltage  remains.  The  boundary  condition  present 
in  BLINE  can  then  be  derived  from  equation  (28b)  of  the  transmission  line 
equations. 
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mm 


3l 

3z 


Gdvd 


(63) 


An  examination  of  figure  3,  which  shows  the  relative  placement  of  the 
current  and  voltage  cells  at  the  source  end  of  the  linear  section  in 
BLINE,  indicates  that  equation  (63)  can  be  differenced  as 


2_  /Tn+1  Tn+1  \ 

Ax  \  j  ‘  1j+h )  'At 


where  the  superscript  designates  the  time  index,  the  subscript  the  position 
index  and 


n+1  _  _1_  Y  n+1  .  hi 
i+h  ~  Dn  X11  Dn 


(65) 


is  the  current-voltage  relationship  at  the  source  end  of  the  linear  section. 
This  gives  the  voltage  at  the  source  end  of  the  linear  section  as 


2_ 

Ax 


hl\ 

Dll/ 


+ 


(66) 


A  similar  analysis  at  the  far  end  of  the  linear  section  shows 


V 


n+1 

j+2 


Gd  2_  !?11  \  „n  (^SL 
T~  "  Ax  Dn  /"  Vj+2\At  '  2 


2_  /.n+1  J_  n+1  J_  n+1  hi  ~n+l\ 
Ax  y  j+2  '  Dn  Xn  ‘  D21  21  '  D21  Vj+l/ 


(67) 


which  uses  a  current  at  the  far  end  of  the  linear  section  of 


I 


n+1 

j+2-% 


N, 


_1_  n+1  +  J1  vn+l  +  _1 
°11  11  D11  3+2  D 


n+1  .  N21  ~n+l  ,,a* 

x2i  +  n-  vi  +  l  (68) 

2!  U2\  31 
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The  voltage  term  V  is  the  voltage  occuring  at  the  source  end  of  the  linear 
•«* 

section.  V  has  been  time  delayed  by  the  length  of  time  it  takes  the  signal 
to  propagate  along  the  linear  section  of  the  cable. 


These  voltage  and  current  relations  were  then  put  into  BLINE  relating 
two  adjacent  voltages.  The  effect  of  this  was  to  separate  the  two  cells 
by  the  length  of  the  inserted  linear  section. 


4.  COMPARISON  WITH  BLINE 

To  compare  the  approximation  of  a  linear  section  by  the  differential 
equations,  the  following  problem  was  run  by  the  new  method  and  BLINE.  A 
coaxial  cable  was  chosen  with  the  parameters 

R  =  .013611  ohms/meter 
L  =  5.8327xl0~8  henries/meter 
G  =  0.0 

C  =  4.3878xl0-10  farads/meter. 

These  were  derived  from  the  physical  parameters 

r,  =  6.35xl0“3  m 

1  -3 

r9  =  8.5x10  m 

c  5 

oc  =  5.8x10  mhos/meter 

e.  =  2.3  farads/meter 
d  o 

where  r^  and  are  the  radius  of  the  core  and  dielectric,  oc  is  the  con¬ 
ductivity  of  the  core,  and  is  the  dielectric  constant  of  the  core  in¬ 
sulation. 

A  transmission  line  of  1.2  km  length  with  these  parameters  was  driven 
at  one  end  with  a  triangular  voltage  pulse  of  one  volt  amplitude  and 
0.5  usee  duration.  For  the  first  comparison,  the  spatial  cell  size  was 
chosen  to  be  10  m  and  observation  points  were  made  at  distances  of  100  m 
and  1100  m.  The  section  of  line  represented  by  the  differential  equations 
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was  1  km  in  length  and  was  inserted  between  two  100  m  sections.  The  ob¬ 
servation  points  were  taken  at  the  two  ends  of  the  1  km  insert.  The  rep¬ 
resentation  of  this  section  used  the  poles  and  residues  discussed  in 
section  III. 3. 

Figure  6  shows  two  voltages:  one  predicted  by  BLINE  100  m  from  the 
source  end  and  the  other  predicted  by  the  modified  BLINE  code  at  the  same 
place.  Figure  7  shows  the  voltages  predicted  by  BLINE  and  the  modified 
BLINE  code  at  the  junction  of  the  insert  1.1  km  from  the  source.  It  can 
be  seen  that  this  agreement  is  not  good  when  BLINE  uses  a  cell  spacing  of 
10  m,  but  improves  when  the  spacing  is  reduced  to  5  m. 

The  discrepancy  is  probably  due  to  the  superior  description  of  the 
transfer  admittance  by  the  differential-equation  technique;  agreement 
probably  occurs  in  the  limit  as  the  cell  spacing  goes  to  zero.  This  is 
shown  by  the  5  m  cell  spacing  resulting  in  a  pulse  more  closely  resembling 
the  pulse  of  the  differential-equation  approach. 
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Figure  6.  Comparison  of  BLINE  and  BLINE  With  ]  km 
Insert  at  the  Source  End  of  the  Insert 


SECTION  V 


LINEAR  SECTION  IN  A  BURIED  CABLE 

To  model  a  realistic  situation  of  cable  breakdown,  two  more  effects 
must  be  included.  The  first  is  to  add  a  voltage  drop  across  the  soil  so 
the  entire  voltage  is  not  across  the  cable  dielectric.  The  second  addi¬ 
tion  is  the  effect  of  a  distributed  source  along  the  linear  insert.  With 
these  additions,  the  regions  of  primary  and  secondary  breakdown  separated 
by  a  linear  length  of  cable  can  be  described  by  the  differential-equation 
technique  and  compared  to  the  predictions  of  the  unmodified  time-domain 
code,  BLINE. 

1.  DISTRIBUTED  SOURCES 

In  the  formulation  described  in  section  II,  the  effect  of  the  dis¬ 
tributed  sources  is  to  add  currents  at  the  ends  of  the  linear  insert 
(equation  (7)).  The  solution  of  the  differential  equations  then  repre¬ 
sents  the  sum  of  the  currents  induced  by  the  distributed  sources  and  the 
currents  resulting  from  signals  propagating  to  the  linear  section 
(equation  (10)). 

In  the  frequency  domain,  it  is  convenient  to  represent  the  E-field 
excitation  of  equation  (2)  as 

e(z,s)  =E  An(s)eian{s)z  (69) 

n 

The  current  sources  will  be  assumed  to  be  zero. 

For  a  general  E-field,  as  a  function  of  time,  the  Fourier  transform 
has  to  be  accomplished  numerically  to  achieve  the  frequency  domain  repre¬ 
sentation  of  equation  (69).  However,  if  the  field  as  a  time  function  is 
a  sum  of  double  exponentials,  the  Fourier  transform  can  be  found  analyti¬ 
cally;  in  this  case, 
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(70) 


E(z,t) 


V  c/9"'2'*0’  (e'a"(t'T") 


-b  (t- 


V 


where  cn  is  the  amplitude,  gn  is  the  spatial  decay  constant,  is  the 
time  delay,  RQ  is  the  data  range,  and  ap  and  are  the  temporal  decay 
constants  of  the  double  exponential.  The  Fourier  transform  of  this 
field  gives 


which  gives 


A„U)  =  cne 


g  R 
an  o 


(71) 


(72a) 


and 


%(s> =  -fe +  % 


(72b) 


by  comparision  with  equation  (69). 

With  the  distributed  source  in  the  form  of  equation  (69),  the  integra¬ 
tion  over  the  length  of  the  linear  insert  can  be  easily  accomplished; 
from  equation  (7),  the  integrals  are 


V1 

Mv1-5'  */ 


T(z  +«-C)  -T(Z  +H-C) 
e  +  e 


E 


lot  C 

A  n  dC 
n 


(73a) 


or 


ia  z 
n  o 


gi<zo’*"s)  * 


n  n 
+  Ysinh(Yi) 


io  l 

-ia  e  +  ia  cosh(Yt) 


(73b) 
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These  show  that  the  first  term  of  Ij,  the  current  addition  at  the  end 
of  the  linear  section  farthest  from  the  source,  has  a  delay  of  ian(zo+£), 
the  light  travel  time  from  the  source  to  the  far  end  of  the  linear  section. 
The  next  term  of  ^  has  a  delay  of  ianzo-m,  the  light  travel  time  to  the 
near  end  of  the  linear  section  plus  the  propagation  time  along  the  linear 
section.  All  other  terms  are  equivalent  to  reflections  at  the  ends  of 
the  linear  sections  with  subsequent  additional  delays  of  2nY£. 

The  first  term  of  I^,  the  current  addition  at  the  end  of  the  linear 
section  closest  to  the  source,  has  a  delay  of  ian20>  the  light  travel  time 
from  the  source  to  that  end  of  the  linear  section.  The  next  term  of  1^ 

has  an  additional  delay  of  (ia  -Y)4;  this  represents  the  delay  of  light 

traveling  from  the  near  to  the  far  end  of  the  linear  section  plus  the 

propagation  delay  along  the  cable  back  to  the  near  end. 

As  in  section  IV,  the  interest  in  predicting  secondary  breakdown  only 
necessitates  the  use  of  those  terms  representing  signals  propagating  away 
from  the  source.  With  this  restriction,  the  sources  can  be  approximated  as 
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2.  ADAPTATION  TO  BLINE 


Two  further  modifications  to  the  code  used  in  the  analysis  of  the 
previous  section  were  necessary.  The  first  involved  the  addition  of  the 
source  terms  1^  and  into  the  boundary  conditions  at  the  ends  of  the 
linear  section.  The  second  was  the  calculation  of  the  dielectric  and 
soil  voltages  at  the  ends  of  the  linear  section.  The  voltage  propagated 
through  the  linear  section  is  the  total  voltage,  the  sum  of  the  dielectric 
and  soil  voltages,  which  necessitates  the  separation  of  the  two  at  the 
boundaries  of  the  linear  and  nonlinear  sections. 


Because  of  time  limitations,  a  case  was  chosen  where  the  soil  con¬ 
ductivity  was  zero;  the  ultimate  reality  of  this  case  was  accomplished 
by  choosing  lumped  parameters  R,  L,  and  C  identical  to  those  of  a  realis¬ 
tic  case  with  frequency  dependent  parameters  at  a  given  frequency;  the 
parameters  were  then  assumed  to  be  constant.  The  reason  for  choosing 
the  soil  conductance  equal  to  zero  was  to  allow  a  simple  relationship 
between  the  dielectric  and  soil  voltages.  From  figure  2,  it  can  be 
seen  that  this  relationship  is 


and 


(78a) 


C  +  C. 
g  d 


(78b) 
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where  V  is  the  total  voltage  and  the  d  and  g  subscripts  refer  to  the 
dielectric  and  ground. 

The  capacitance  across  the  ground,  C^,  was  calculated  by  assuming  that 
the  total  capacitance,  C,  of  the  realistic  calculation  was  the  result  of 
the  soil  and  dielectric  capacitances  in  series. 

The  boundary  conditions  at  the  junction  of  the  linear  and  nonlinear 
sections  are  modifications  of  equations  (66)  and  (67)  including  the  cur¬ 
rent  source  and  voltage  fraction  V^/V .  The  dielectric  and  ground  voltages 
were  found  from 


at  the  end  of  the  linear  section  nearest  the  source  and 
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at  the  other  end  of  the  linear  section.  In  both  cases 
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The  source  terms,  S1  and  Sj,  are  the  inverse  Fourier  transforms  of 
equations  (77a)  and  (77b).  All  other  terms  are  identical  to  those  of 
equations  (66)  and  (67). 
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3.  COMPARISON  WITH  BLINE 

The  comparison  problem  to  be  done  by  BLINE  and  the  differential  equa¬ 
tion  method  was  determined  by  first  finding  the  lumped  parameters  of  a 

5 

buried  cable  at  a  frequency  of  1.45  x  10  Hz.  The  frequency  chosen  was 
an  estimate  of  the  frequency  of  the  pulse  which  is  responsible  for  the 
generation  of  secondary  breakdown.  Even  though  the  conductance  was  non¬ 
zero  at  this  frequency,  its  effects  were  ignored.  The  input  parameters 
were 

oc  =  5.8  x  107  mhos/meter 

r^  =  6.35  x  10"^  m 

r^  =  8.5  x  IQ-"*  m 

=  2.3  farads/m 

%  water  content  of  the  soil  3  10%. 

5 

At  a  frequency  of  1.45  x  10  Hz,  the  lumped  transmission  line  parameters 

were 

R  =  1.4878  x  10" 1  ohms/m 
L  =  1.4906  x  10"®  henries/m 
C  =  4.3379  x  10"10  farads/m. 

These  parameters  were  then  used  as  constant  R,  L,  and  C  for  the  linear 
section  of  the  line  to  be  represented  by  the  differential-equation  tech¬ 
nique. 

In  the  sections  on  either  end  of  this  linear  region  and  the  entire 
cable  represented  by  BLINE,  the  voltage  is  divided  into  voltage  across 
the  dielectric  and  voltage  across  the  ground;  in  this  case  the  parameters 
are 


* 
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where 


R  =  1.4878  x 
L  =  1.4906  x 
Cd  =  4.3878  x 
Cg  -  3.8175  x 


10" 1  ohms/m 
10~8  henries/m 
10"10  farads/m 
10"8  farads/m 


Cd  =  2lT  ed  /  ln  (ri/r2^ 


and  Cg  is  calculated  from 


Using  these  parameters  and  a  dielectric  breakdown  voltage  of  107  v/m, 
BLINE  was  run  to  find  the  breakdown  regions  of  a  4  km  cable.  The  network 
was  terminated  at  the  source  end  by  a  10  ohm  resistor  and  at  the  far  end 

4 

by  a  10  ohm  resistor;  this  mismatch  was  purposely  chosen  to  enhance 
secondary  breakdown. 

With  this  input,  BLINE  predicted  a  primary  breakdown  region  of  2.6  km 
and  a  secondary  breakdown  region  of  60  m,  4  km  from  the  burst. 

The  region  from  2.8  to  3.8  km  remained  linear.  This  was  chosen  as  the 
section  to  be  replaced  by  a  line  represented  by  the  differential -equation 
technique. 

To  obtain  the  differential  equations,  the  Prony  analysis  was  done 
on  the  inverse  Fourier  transforms  of  the  self  and  transfer  admittances. 

As  in  the  case  discussed  in  section  IV,  the  admittances  were  limited  to 
the  terms  representing  the  signals  traveling  away  from  the  source.  The 
poles,  residues,  and  impulse  terms  for  the  self  admittance  were 
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aj  =  1.63148  x  102  a2  =  6.88527  x  102 

sl  =  -2.5697  x  104  s2  =  -8.64688  x  104 

Impulse  =  1.7059  x  10”2. 

The  terms  representing  the  transfer  admittance  were 

aj  =  1.14429  x  102  a2  =  6.04022  x  101 

s,  =  -8.20423  x  104  s2  =  -1.421511  x  104 

Impulse  =  9.5908  x  10"3  . 

As  in  the  previous  case  (section  IV. 2),  a  satisfactory  fit  was  achieved 
for  both  admittances  with  two  poles  and  residues.  Because  the  breakdown 
phenomenon  involves  longer  cables,  the  Prony  fits  were  found  to  a  time  of 
60  ys. 

The  computation  of  the  source  terms  was  done  in  the  frequency  domain 
and  inverse  Fourier  transformed  at  times  that  could  be  directly  utilized 
by  the  time-domain  code  at  the  two  junction  points. 

To  compare  the  two  processes,  observations  were  made  of  the  dielectric 
voltage  at  the  two  junction  points,  2.8  km  and  3.8  km  along  the  cable,  and 
at  the  end  of  the  line  where  secondary  breakdown  was  expected.  Because 
the  primary  breakdown  region  is  essentially  described  by  the  same  coding, 
no  differences  in  this  region  were  found  between  the  two  predictions. 

Figure  8  shows  the  responses  predicted  at  the  junction  of  the  primary 
breakdown  region  with  the  linear  section  of  cable.  The  dotted  line  (at 
70  ys),  in  this  figure  and  the  following  two,  is  the  limit  to  which  the 
Prony  approximations  of  the  admittances  were  made.  At  larger  times,  an 
extrapolation  is  indicated.  In  figure  8,  the  variation  between  the  pre¬ 
dicted  responses  at  80  ysec,  is  the  result  of  a  reflection  from  the  end 
of  the  cable;  because  only  signals  traveling  from  the  source  were  handled 
by  the  differential  equations,  the  reflection  is  not  seen  with  this  analy¬ 
sis.  Figure  9  compares  the  response  predicted  by  BLINE  3.8  km  from  the 
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Across  the  Dielectric  Voltage  Across  the  Dielectric 


source  with  the  response  generated  by  the  differential  equations.  The 
last  comparison  between  the  two  approaches  is  shown  in  figure  10.  This 
shows  the  predicted  responses  at  the  end  of  the  cable  4  km  from  the  source 
The  oscillations  are  due  to  reflections  off  the  end  of  the  cable  and  the 
secondary  breakdown  region. 


24x10 


Time  (ys) 


Figure  10.  Dielectric  Voltages  Predicted  by  BLINE  and  the  Differential 
Equation  Technique  at  the  Cable  Termination 
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The  starting  time  and  extent  of  secondary  breakdown  predicted  by  the 
two  methods  are 

Time  Extent 

BLINE  49.54  ps  40  m 

Diff.  Eqn.  49.34  ps  60  m 

The  earlier  time  and  greater  extent  of  secondary  breakdown  predicted 
by  the  differential  equation  approach  is  probably  a  result  of  a  more 
accurate  representation  of  the  propagation  of  the  pulse  through  the  1  km 
linear  region.  A  similar  difference  in  pulse  propagation  was  shown  in 
the  coaxial  line  comparison  in  section  IV. 4. 
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SECTION  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 

1.  CONCLUSIONS 

Two  methods  have  been  employed  to  represent  a  linear  section  of  cable 
between  the  regions  capable  of  describing  nonlinear  behavior. 

The  two  methods  employ  a  relationship  between  the  currents  and  voltages 
at  the  ends  of  the  linear  section.  The  first  uses  convolution  integrals 
to  relate  the  frequency-domain  response  of  current  and  voltage  and  presents 
the  relationship  in  the  time  domain.  The  disadvantage  of  this  approach 
is  the  necessity  of  the  calculation  of  convolution  integrals  at  each  time 
step.  The  second  approach  recognizes  that  the  convolution  integrals  in 
this  case  are  the  solutions  of  differential  equations.  Once  the  differential 
equations  are  obtained,  each  time  step  requires  the  addition  of  one  incre¬ 
ment  to  the  differential -equation  solution  rather  than  the  calculation  of 
an  entire  convolution  integral. 

The  use  of  convolution  integrals  has  shown  the  ability  of  this  method 
to  describe  the  self  and  transfer  admittances  of  a  lossless  line  and  account 
for  the  reflections  of  pulses  in  the  network. 

The  differential  equation  method  has  satisfactorily  described  the  prop¬ 
agation  of  a  pulse  through  a  1  km  section  of  coaxial  cable.  In  all  prob¬ 
ability,  the  pulse  exiting  the  cable  described  in  this  manner  is  closer  to 
reality  than  the  pulse  described  by  a  pure  time-domain  approach. 

This  procedure  has  also  been  used  to  predict  the  secondary  breakdown 
occurring  in  a  buried  cable  4  km  in  length.  The  time  and  extent  of  the  pre¬ 
dicted  nonlinearities  and  the  breakdown  described  by  BLINE  are  essentially 
in  agreement. 
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2.  RECOMMENDATIONS 


The  use  of  constant  R,  L,  G,  and  C  parameters  to  describe  the  linear 
section  is  most  likely  inadequate  for  most  applications.  However,  the 
use  of  frequency  dependent  parameters  to  describe  the  admittances  is  dif¬ 
ficult  because  of  the  lack  of  well  defined  impulse  behavior.  It  is  rec¬ 
ommended  that  further  study  be  made  of  the  three  media  problem  to  determine 
the  high  frequency  limits  of  the  self  and  transfer  admittances  of  a  linear 
transmission  line. 

The  technique  of  using  differential  equations  to  describe  convolution 
integrals  can  be  applied  to  each  spatial  cell  of  a  nonlinear  time  domain 
code  to  achieve  a  significant  increase  in  the  predictive  capability  of  such 
a  code.  This  coupled  with  the  previous  recommendation  would  enable  the 
time-domain  code  to  be  as  accurate  as  a  frequency-domain  code  but  still 
enable  it  to  describe  nonlinearities. 

An  easily  achieved  improvement  on  the  present  use  of  differential 
equations  in  this  report  is  the  addition  of  soil  conductivity  to  the  pro¬ 
cedure.  A  simple  relationship  between  the  total  voltage  propagated  by 
the  linear  section  and  the  dielectric  voltage  needed  at  the  junction  to 
the  nonlinear  section  is 

3V. 

< W  vd +  <Vcg>  it 

where  the  d  and  g  subscripts  refer  to  the  dielectric  and  ground  respectively 
and  V  is  the  total  voltage,  V^+Vg. 
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